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Abstract. We study the influence of convection on the asymmetries and Doppler shifts 
of Fei spectral lines in the metal-poor red giant HD 122563. To this end, we compute the- 
oretical Fei line shifts and line bisectors using 3D hydrodynamical model atmosphere of 
HD 122563 calculated with the C0 5 B0LD code. We then make a detailed comparison of the 
theoretical line shifts and bisectors with those derived from the high quality HARPS spec- 
trum of HD 122563 taken from the ESO Science Archive Facility (R = 115 000, average 
signal-to-noise ratio, S/N ~ 310). In general, we find a good agreement between the the- 
oretically predicted and observed Doppler shifts of Fe i line cores, with somewhat larger 
discrepancies seen in the case of weaker (equivalent width W < 5 pm) and stronger lines 
fT") ' (W > 1 1 pm). Both observed and theoretical coreshifts cover a range between and -1 km/s, 

with increasingly stronger blueshifts for weaker lines and slight hints of a coreshift depen- 
\Q , dence on wavelength. Theoretical bisectors reproduce the observed ones reasonably well 

too, however, theoretical bisectors of the weak red (A > 600 nm) Fei lines have blueshifts 
that are by up to ~ 200 m/s larger than observed. The obtained results therefore suggest that 
the current C0 5 B0LD models are capable of reproducing the large-scale velocity fields in 
the atmosphere of HD 122563 sufficiently well. Nevertheless, further efforts are needed in 
order to understand the physical reasons behind the discrepancies in theoretical predictions 
and observed properties of the weakest and strongest Fe i lines. 
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Convection is an important physical phe- 

nomenon responsible for the energy transport 
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stars (such as the Sun), convection zone may 
reach close to (or even into) the atmospheric 
layers where the formation of spectral lines 
takes place, and therefore may directly influ- 
ence line strengths, shapes, and Doppler shifts. 
A direct imprint of convection is the granu- 
lation pattern seen on the surface of the Sun, 
comprised of hot rising granules and cooler in- 
tergranular downflows. Granulation cannot be 
directly resolved in other stars yet, with the 
exception of the closest supergiants such as 
Betelgeuse (e.g., Kervella et al. 2009), there- 
fore it can be only detected indirectly, for ex- 
ample, through the line asymmetries and their 
Doppler shifts. 

3D hydrodynamical stellar model atmo- 
spheres, such as those calculated with the 
Stagger (Stein & Nordlund 1998) and 
C0 5 B0LD (Freytag et al. 2012) codes, treat 
convection in considerably more realistic way 
than the classical ID models do. Up to now, 
Stagger and C0 5 B0LD models were success- 
fully applied to study spectral line properties 
in the Sun (Asplund et al. 2000; Caffau et al. 

2011) , Procyon (Allende Prieto et al. 2002), a 
main sequence turn-off star HD 74000 (Cayrel 
et al. 2007), several M dwarfs (Ramirez et 
al. 2009), and red giants (Collet et al. 2007; 
Ramirez et al. 2010; Kucinskas et al. 2013; 
Dobrovolskas et al. 2013). These studies have 
shown that theoretical predictions were reason- 
ably consistent with the observations, suggest- 
ing that current 3D hydrodynamical model at- 
mospheres are already sufficiently realistic to 
account for the convective flows taking place 
in real stellar atmospheres. 

In this work we present first results of our 
further application of the 3D hydrodynamical 
C0 5 B0LD model atmospheres to study convec- 
tion in red giant stars. We focus our investiga- 
tion on the asymmetries and Doppler shifts of 
Fe i lines in the atmosphere of the metal-poor 
giant HD 122563. Atmospheric parameters of 
this star have been re-derived recently by us- 
ing HIPPARCOS parallaxes to obtain its surface 
gravity, log g = 1 .6 (van Leeuwen 2007), and 
CHARA interferometry to determine its effec- 
tive temperature, T e g = 4600 K (Creevey et al. 

2012) . Spectral line coreshifts and bisectors in 
HD 122563 have been studied recently with the 



aid of Stagger models (Ramirez et al. 2010). 
All this makes HD 122563 an interesting tar- 
get for a study with CO 5 BOLD models, allowing 
to test the realism of the C0 5 B0LD model atmo- 
spheres in the metal-poor giant regime, and to 
compare theoretical predictions obtained with 
the Stagger and C0 5 B0LD codes. 

2. Observed spectrum of HD 122563 

In this study we used a publicly available 
reduced HARPS spectrum of HD 122563 
extracted from the ESO Science Archive 
Facility 1 (program ID 080.D-0347(A)). The 
spectrum was obtained in the wavelength range 
of 380 - 680 nm, with the average signal-to- 
noise ratio of S/N « 310 and spectral reso- 
lution of R = 115000. We used the pipeline- 
reduced spectrum which was processed with 
the HARPS Data Reduction Software 2 (DRS). 
The continuum normalization was done using a 
second order polynomial fit with the Dech28T 
spectral analysis package 3 . 

3. Three-dimensional C0 5 B0LD model 
atmospheres and spectral line 
synthesis 

For the calculation of synthetic Fe i lines we 
used two C0 5 B0LD model atmospheres calcu- 
lated with identical effective temperature and 
gravity (corresponding to those of HD 122563, 
T eS = 4590 K and logg = 1.6, cgs), at two dif- 
ferent metallicities, [M/H] = -2.0 and -3.0. 
An enhancement of 0.4 dex was applied to the 
abundances of O, Ne, Mg, Si, S, Ar, Ca, and 
Ti, while for other elements solar-scaled abun- 
dances from Grevesse & Sauval (1998) were 
used (solar-scaled abundances from Asplund et 
al. 2005, were used in case of C, N, and O). The 
C0 5 B0LD model atmospheres were calculated 
using a rectangular box of 200 x 200 x 170 grid 
points and covering 4.2x4.2x2.1 Gm in phys- 
ical dimensions (x, y, z, respectively). From the 

1 http : // archive . eso . org/eso/eso_ 
archive_main . html 

2 http : // archive . eso . org/ archive/adp/ 
ADP/HARPS/ index . html 

3 http://www.gazinur.com/ 
DECH- software . html 
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entire 3D model run, twenty 3D model snap- 
shots (i.e., 3D model structures obtained at dif- 
ferent instants in time) were chosen for the 
spectral synthesis of the Fei lines. The snap- 
shots were selected in such a way as to en- 
sure that properties of the 20 snapshot ensem- 
ble would match those of the entire 3D model 
run as closely as possible. 

Synthesis of the Fei lines was performed 
under the assumption of local thermodynamic 
equilibrium (LTE), using for this purpose 
the Lin£or3D spectral synthesis package 4 . 
Spectral lines were synthesized for different 
iron abundances to obtain sequences of Fei 
line profiles at [M/H] = -2.0 and -3.0, with 
the iron abundances in steps of 0.2 dex. The ob- 
tained line profiles were then interpolated for 
the metallicity of HD 122563, [Fe/H] = -2.6 
(the latter taken from Gratton et al. 1996). To 
speed up the line synthesis calculations, only 
every third point of the 3D model structure (in 
each horizontal direction) was used to synthe- 
size the Fei lines. Tests made with the full 3D 
model structures showed that this simplifica- 
tion does not alter the resulting line bisectors 
by more than ~ 0.02 km/s. 

For the line bisector and coreshift analy- 
sis (Sect. 4) we used 8 1 non-blended Fe i lines 
located in the wavelength range of ~ 400- 
670 nm. Line wavelengths were taken from 
Nave et al. (1994), the excitation potentials 
and line broadening parameters were extracted 
from the VALD database (Kupka et al. 1999). 
The rotational velocity (vsini) was obtained 
by simultaneously fitting the line strength, 
Doppler shift, and v sin i of the theoretical and 
observed line profiles, keeping the instrumen- 
tal broadening and continuum level fixed. 



4. Line asymmetry measurements 

Three spectral line properties were used to 
study the influence of convection on the for- 
mation of Fe i lines: equivalent width (W, used 
as line strength indicator), line coreshifts (i.e., 
Doppler shifts of the deepest points in the line 



4 http : //www. aip . de/~mst/Linfor3D/ 
linfor_3D_manual .pdf 



profiles), and line bisectors 5 . Equivalent widths 
were measured by integrating theoretical line 
profiles obtained as best-fits to the given ob- 
served line profile. Both observed and theoreti- 
cal coreshifts and bisectors were calculated us- 
ing the following identical routines: 

- coreshifts were computed by fitting 
Gaussian profile to the five points cen- 
tered around the deepest point in the 
observed/synthetic spectral line profile; 

- bisectors were calculated utilizing the 
usual procedure of finding midpoints be- 
tween the red and blue line wings at a 
number of different line depression lev- 
els. To obtain bisector values at any given 
flux value, a monotonized interpolation 
was used to interpolate between the ob- 
served/synthetic points on the red and blue 
wings of the line profile. 

Errors in the observed line coreshifts were 
computed by adding uncertainties due to lab- 
oratory wavelengths, HARPS wavelength cal- 
ibration, and photon noise (converted into 
wavelength uncertainty) in quadrature and tak- 
ing a square root of the sum. We assumed that 
the error in the laboratory wavelength calibra- 
tion was 2m/s (taken from Nave et al. 1994) 
and the error in the wavelength calibration « 
15 m/s (HARPS DRS estimate). The contribu- 
tion from the photon noise to the coreshift error 
was estimated using the x 2 fitting procedure. 
For this, we assigned the measurement error to 
each flux point in the fitted line profile (cal- 
culated from the local S /N), fitted it with the 
Gaussian profile, and obtained the coreshift er- 
ror corresponding to the x 1 ± 1 brackets. One 
should note that uncertainty due to the pho- 
ton noise depends both on the line strength and 
wavelength (S /N increases towards the red part 
of the spectrum). 

In case of the observed line bisectors, we 
accounted only for errors due to photon noise 
which was computed using the formula of 
Gray (1983). Uncertainties due to the labo- 
ratory wavelength measurements and wave- 

5 The bisector is a line connecting loci of the mid- 
points between the blue and red wings of a spectral 
line profile at different flux (line depression) levels. 
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Fig. 1. Top panel: observed and theoretical coreshifts of Fe i lines in the metal-poor red giant 
HD 122563, plotted versus line equivalent width, W. Errors of the observed coreshifts take into 
account uncertainties in the laboratory wavelength, instrumental wavelength calibration, and 
photon noise. Error bars of the theoretical coreshifts show l<x of the coreshift distribution ob- 
tained from the Monte Carlo bootstrap experiment (see text for details). Bottom panel: difference 
between the observed and theoretical coreshifts (i.e., residual coreshifts) plotted versus W. 



length calibration errors were not taken into ac- 
count since they do not change across the given 
spectral line profile. 

It is important to remind that theoretical 
line profiles obtained using each of the twenty 
3D model snapshot selection are in fact dif- 
ferent from each other. Moreover, the proper- 
ties of spectral lines computed using any dif- 
ferent ensemble of the twenty 3D model snap- 
shots should be slightly different from those 
obtained using the current snapshot selection. 
We therefore applied a bootstrap Monte Carlo 
method to estimate statistical uncertainty in the 
bisectors and coreshifts of the synthetic spec- 
tral lines. To this end, from the initial pool of 
twenty spectral lines computed using twenty 
3D model snapshots we randomly selected 
twenty line profiles (allowing each profile to be 
selected several times). This procedure was re- 
peated 1000 times, resulting in 1000 different 
sets of twenty theoretical line profiles. We then 
computed average core shifts and bisectors for 



each of these 1000 sets of spectral line pro- 
files, and calculated their standard deviations 
(in the case of bisectors we computed standard 
deviations at a number of different flux lev- 
els). These standard deviations were then used 
as uncertainties in the theoretical line bisectors 
and coreshifts. 

5. Results and discussion 

Before comparing observed and theoretical 
coreshifts, we determined the radial velocity 
(vrad) of HD 122563 by subtracting theoret- 
ical coreshifts from the observed ones. The 
mean residual was then added to the mean ob- 
served line shift to obtain the radial velocity, 
Vrad = -25.57 ±0.11 km/s (the error represents 
the standard deviation). An identical approach 
was utilized by Ramirez et al. (2010) who ob- 
tained a RV v rad = -25.39 + 0.09 km/s us- 
ing 3D hydrodynamical Stagger atmosphere 
models. The agreement between the two val- 




Fig. 2. Top panels: examples of observed and best fitting theoretical Fe i line profiles in the metal- 
poor red giant HD 122563. Bottom panels: bisectors of the spectral lines shown in the top panel. 
Gray shade represents the observed bisector uncertainty (lcr), red shaded region shows the lcr 
uncertainty of the theoretical bisector. 



ues obtained with different 3D hydrodynamical 
model atmospheres is indeed very good (the 
difference is 0.16 km/s). 

The observed coreshifts (corrected for the 
radial velocity of HD 122563) generally agree 
well with those predicted using the CO 5 BOLD 
models (Fig. 1). There are hints for the core- 
shift dependence on wavelength seen both in 
the observed and theoretical coreshifts, and 
these trends may be partly responsible for the 
scatter in the W-coreshift plane seen at any 
given value of W. The agreement between the 
observed and theoretical coreshifts is reason- 
ably good in the case of medium strong lines 
(5 pm < W < 11 pm), with the mean resid- 
ual coreshift of 0.015 + 0.076 km/s (error is 
lcr scatter around the mean value). The situ- 
ation is slightly worse in the case of weakest 
(W < 5 pm) and strongest (W > 1 1 pm) lines, 
with the mean residual coreshifts of +0.09 + 
0.13km/s and -0.088 + 0.081 km/s, respec- 
tively. These discrepancies are about two times 
larger than the respective median error of the 



individual residual coreshift, 0.044 km/s (weak 
lines) and 0.062 km/s (strong lines). 

Spectral line bisectors may be yet another 
tool to study velocity fields in stellar atmo- 
spheres, especially since they carry the im- 
prints from the velocity fields and tempera- 
ture fluctuations over the formation region of 
a given spectral line. In our study, observed 
and theoretical bisectors were compared by 
correcting theoretical bisectors for the radial 
velocity of HD 122563. We stress, however, 
that no additional corrections were applied to 
the theoretical bisectors, as it is sometimes 
done in other bisector studies (e.g., Ramirez 
et al. 2010). One may note that in case of 
both observed and theoretical bisectors there 
is a clear relation between the strength of a 
given spectral line and its bisector shape, with 
the bisectors of strong lines being nearly verti- 
cal and those of weaker lines becoming more 
curved, with the line cores getting progres- 
sively more blueshifted with decreasing line 
strength (Fig. 2). 
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From the comparison of observed and the- 
oretical bisectors provided in Fig. 2 we con- 
clude that there is a reasonably good agree- 
ment between the observed and theoretical bi- 
sectors of strong lines and lines of intermediate 
strength (apart from the residual discrepancy 
in the coreshifts of the strong lines). However, 
in the case of weak lines discrepancies in both 
the coreshifts and bisector shapes are more 
pronounced. This result is similar to that ob- 
tained by Ramirez et al. (2010) with the aid 
of Stagger models, where the weakest theo- 
retical bisector was blueshifted by 0. 1 km/s to 
match the observed bisector (the blueshifts re- 
quired in case of stronger lines were smaller, 
at the level of ~ 0.03 km/s). Reasonably good 
agreement between theoretical predictions and 
observations may therefore suggest that cur- 
rently available C0 5 B0LD models are reproduc- 
ing the global velocity fields in the atmosphere 
of HD 122563 sufficiently well. Nevertheless, 
further improvements may be needed in order 
to obtain better agreement in case of the weak- 
est and strongest spectral lines. 



6. Conclusions 

Preliminary results obtained during the project 
aimed to study spectral line asymmetries in the 
metal -poor giant HD 122563 show that Fei 
line coreshifts and bisectors predicted by the 
current 3D hydrodynamical C0 5 B0LD model at- 
mospheres agree reasonably well with those 
observed in HD 122563, despite the minor 
problems in case of the weakest and strongest 
spectral lines. These results are compatible 
with those obtained in a similar study of line 
asymmetries in HD 122563 carried out us- 
ing 3D hydrodynamical Stagger model atmo- 
spheres, pointing to the consistency in the pre- 
dictions made with two different 3D hydrody- 
namical model atmosphere codes. All this sug- 
gest that the current level of realism of the 3D 
hydrodynamical model atmospheres should be 
sufficient to warrant further applications of the 
3D model atmospheres in a wider variety of 
astrophysical contexts, in order to further test 
their realism and identify problems where ad- 
ditional improvements may be needed. 
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